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The time dependent density matrix of a system with potential barrier is studied using path 
integrals. The characterization of the initial state, which is assumed to be restricted to one side of 
the barrier, and the time evolution of the density matrix lead to a three-fold path integral which is 
evaluated in the semiclassical limit. The semiclassical trajectories are found to move in the complex 
coordinate plane and barrier penetration only arises due to fluctuations. Both the form of the 
semiclassical paths and the relevant fluctuations change significantly as a function of temperature. 
The semiclassical analysis leads to a detailed picture of barrier penetration in the real time domain 
and the changeover from thermal activation to quantum tunneling. Deep tunneling is associated 
with quasi-zero modes in the fluctuation spectrum about the semiclassical orbits in the long time 
limit. The connection between this real time description of tunneling and the standard imaginary 
time instanton approach is established. Specific results are given for a double well potential and an 
Eckart barrier. 
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I. INTRODUCTION 



Semiclassical theories have been found extremely powerful in understanding the dynamics of complex quantum 
mechanical systems. Special attention has been paid to theories of tunneling processes as they occur in physics, 
chemistry, and biology. Currently, a variety of quantum rate theories are in use explaining experimental findings 
for several situations of interest [GJ. Thereby, roughly speaking, two different strategies can be distinguished. The 
first class of approaches constructs the rate from purely thermodynamic considerations. An example is the bounce 
or instanton method (also called ImF method), originated by Langer Q and extended by several authors B. In 
essence, tunneling rates are derived from the imaginary time dynamics in the inverted potential. Other approaches 
of this type |4|5|] start from periodic orbit theory ](| in imaginary time. Tractable rate formulae are obtained with 
the centroid method j7j] leading to a semi-empirical separation of dynamical and thermal factors. These methods are 
computationally very efficient and have been applied successfully to systems as diverse as tunneling centers in metals, 
Josephson junctions, or hydrogen bonds, to name but a few. However, closer examination reveals that these theories 
are based in one way or another on ad hoc assumptions that are not derived from first principles. For instance, the 
ImF method postulates a relation between the decay rate and the imaginary part of the free energy. In fact, in some 
cases thermodynamic methods fail to predict the correct rate, e.g. they do not reproduce the energy diffusion limited 
decay for very weakly damped systems at finite temperature. Moreover, these methods are designed to describe 
incoherent decay only. 

The second class of theories describes barrier crossing in terms of dynamical quantities. Perhaps most familiar is 
Yamamoto's rate formula [pi, in essence a Kubo type formula relating the rate with a flux- flux correlation function. 
As shown in Q it is exact only for scattering problems, while in multi- or metastable systems one has to assume 
the existence of a plateau region for times long compared to typical relaxation times but small compared to decay 
times. This restricts the approach to incoherent decay. Recent advances in the real time description of barrier crossing 
based on flux-flux correlations were made e.g. by Voth et al. ]To| ], who incorporated analytically known dynamical 
factors into the rate expressions. Quite recently, Pollak and coworkers jll],[l2| were able to improve this idea by using 
a thermally symmetrized flux operator. A different way of including dynamical information in an approximate way, 
favored by Miller and coworkers, employs semiclassical Initial Value Representation (IVR) for the quantum propagator 
jb|. Although quite successful at high to moderate temperatures, where the quantum dynamics is governed by quasi- 
classical above-barrier processes, these approaches usually fail at low temperatures where deep tunneling prevails. 
Finally, we mention a kind of hybrid approach, the "real time" instanton theory JLlj which includes tunneling in 
the real time propagator by means of instantaneous tunneling transitions. This method is restricted to multi-stable 
systems in the low temperature limit ]l5| . Hence, although the semiclassical theory of quantum tunneling is often 
regarded to be well settled, this turns out to be true only for some limiting cases. What would be desirable is a 
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semiclassical theory starting from first principles that covers the entire range of temperatures as well as coherent and 
incoherent tunneling processes. 

In the realm of classical physics the theory of thermally activated rates is rather firmly based. In a seminal paper 
p6| Kramers determined thermal decay rates from the equation of motion for the phase space distribution function, 
i.e., from the real time dynamics of the system. A corresponding treatment of tunneling in the semiclassical limit 
seems not to be possible, since all real time minimal action paths connecting two sides of the barrier have energies 
larger than the barrier energy. In the dynamical approaches discussed above these trajectories do account for tunneling 
corrections to classical rates [ [l3|]l7| ] , but it is usually argued Jl^,|ll| that within a semiclassical theory deep tunneling 
can only be described by incorporating in addition imaginary time trajectories as they are used in thermodynamic 
methods. The lack of a first principles semiclassical theory of tunneling is intimately connected with a fundamental 
shortcoming of semiclassical real time propagators derived from a single dominant path. Since tunneling arises from 
coherent interference of waves, a satisfactory semiclassical theory needs to capture this interference pattern in terms of 
an appropriate family of real time paths rather than making ad-hoc modifications of simple semiclassical propagators. 

Very recently, we have proposed a theory for transport across a barrier based solely on the real time dynamics 
of the density matrix and the semiclassical approximation |2(J. As a notable feature, the method applies equally 
well to damped and undamped systems and comprises in a unified way the entire range from thermally activated 
decay to low temperature tunneling. In appropriate limits the results of other methods are recovered. Here, we 
explain the technicalities of the approach and evaluate the semiclassical propagator in detail for two paradigmatic 
systems, namely, a double well potential and an Eckart barrier. Thereby, we mainly concentrate on one-dimensional 
models. It turns out that the corresponding theory already reveals the basic structure and that the generalization to 
multidimensional systems, though tedious in detail, is straightforward within the path integral formalism [pl| . 

The article is organized as follows. Next, in Sec. 0, we outline the general semiclassical theory which is then used 
in Sec. [II to derive as a simple example the stationary flux across a parabolic barrier. The main part of the paper 
studies the real time dynamics from high down to vanishing temperature for the cases of a double well potential 
(Sec. IV) and an Eckart barrier (Sec. [v]). Finally, in Sec. VI we summarize the main features of the approach and 
present our conclusions. 



II. GENERAL THEORY 



We consider a statistical ensemble of quantum mechanical particles of mass M moving in a barrier potential V(q) 
at inverse temperature (5 — \/k^T. We chose the coordinate q so that the barrier top is located at q = and measure 
energies relative to the barrier energy by putting V(0) = 0. The initial nonequilibrium state is assumed to be of 
the form of an equilibrium state restricted to the left side of the barrier. Below we will invoke the semiclassical 
approximation which is appropriate provided the barrier height V& is by far the largest energy scale in the system. 
The time evolution of the density matrix p(t) = exjp(-iHt/h)p(0)exp(iHt/h) reads in coordinate representation 

p{lf,Qf,t) = j dqidq' i Gt{qf,q i )p{qi,q' i ,0)Gt(q' f ,q^)* (1) 

where the real-time propagator is given by 

G t (q,q')=(q\e X p(-iHt/h)\q') (2) 

and p(9i,9i,0) describes the initial state. In principle, for our purpose any initial distribution that matches onto 
equilibrium on the left side and vanishes on the right side of the barrier top is appropriate. As long as the restricted 
equilibrium state gives vanishing probability to find the particle on the right side of the barrier top, different initial 
preparations lead to the same long time behavior of the density matrix. Here we put explicitly 

p(q, q', 0) = Z-'ppiq, q>)d(-q)8(- q >) (3) 

for convenience with the proper normalization factor Z and the equilibrium density matrix 

P0 {q,q') = {q\eM-PH)\q'). (4) 

Now, employing the path integral representation for exp(±itH/h) and exp(— /3-ff), respectively, the above integrand 
in Eq. (Q) can be written as a three- fold path integral where two real time paths q(u) and q'(u) run in the interval 
< u < t from qi and q[ to fixed endpoints g/ and q'j, respectively, while those former coordinates are connected by 



2 



an imaginary time path q(a) in the interval < a < h(3, see fig. 0. The real time paths describe the time evolution of 
the system and the imaginary time path the initial state. Of course, a more complete theory would explicitly include 
the coupling to a heat bath environment. In fact, the general scheme of this approach in the case of damped systems 
has already been given elsewhere pT| . Much of the analysis presented below can in principle be extended to this 
situation, however, only a limited number of steps can be carried out analytically due to the more complicated form 
of the effective action functionals. Here, we limit ourselves to undamped motion which allows to treat deep tunneling 
without resorting to numerical methods. This way, the guiding concepts will become more transparent. 

The density matrix p(q,q',t) contains all information about the noncquilibrium quantum process, in particular, 
the average of the operator F — [pS(q) + 5{q)p]/2M gives the flux out of the metastable state, i.e., in coordinate 
representation 



J(t) = (h/2iM) [dp(q f ,~q f ,t)/dq f } qf=0 . 



(5) 



If the flux becomes quasi-stationary, J(t) = Jfl within a certain "plateau region" of time, the escape rate follows from 

r = j fl . 

While an exact solution of Eq. (]]]) for anharmonic barrier potentials is not possible, a high barrier naturally suggests 
a semiclassical approximation. In the semiclassical limit the above path integrals are dominated by minimal action 
paths determined by Hamilton's equation of motion either in the potential V(q) (for the real time propagators) or 
—V(q) (for the equilibrium density matrix). Each path contributes with an exponential factor containing its minimal 
action and a prefactor arising from the Gaussian fluctuations about the minimal action paths. Specifically, the action 
in real time reads 



S(q,q') 



du[Mq 2 /2 - V(q)} 



(0) 



while its imaginary time version, the so-called Euclidian action, is given by 

S(q,q')= / da[Mq 2 /2 + V(q)}. 
Jo 



(7) 



Thus, in the Gaussian semiclassics the propagator (121) is approximated as 



G t (q,q')= V^TIexp 
cl. paths 



ls{q, q ')-i\v 



(8) 



where A(q, q 1 ) = [— d 2 S(q, q') / dqdq']/2iTih and v is the Maslov index. Throughout this paper we also use an equivalent 
representation of the prefactor, namely, 



A(q,q') 



iM 
2nh 



g(°)gW Q E 2 



(9) 



where W(q,q') = f dq"p — S(q,q') + Et is the short action. The corresponding approximation to the equilibrium 

density matrix (||) follows by formal analytic continuation t — ► —ih/3, i.e., by replacing S(q,q') by iS(q,q') in Eq. (||) 
with v = 0. As a result, the integrand in Eq. ([!]) is completely determined by classical mechanics in real and imaginary 
time, respectively, and dominated by an action factor exp[— S(gj, q'f\qi, q'i)/Ti — in(v — v')/2] with 



s (9/,9/l*,<70 = -iS{q f ,q t ) + S(q t ,q' t ) +«S , (gj,<^). 



(10) 



With the approximate integrand at hand, it is consistent to evaluate the ordinary integrations in Eq. (Q) in stationary 
phase. The stationary phase points are determined by minimizing S with respect to the initial coordinates qi-,q'i, i.e., 



<9S 



= 0, 



(?/, q' f ) 



9E 



= o. 



(is-i' f ) 



(11) 



Since the endpoints qf,q'f are fixed, the resulting stationary phase points q s (t) and q' s (t) are functions of time with 
q s (0) — qf, q' s (0) = q'f. For finite t these roots are in general complex. The dominant imaginary time path q s (a) 
connects q' s (t) with q s (t), and the two real time paths q{u) and q'(u) connect q s (t) and q' s (t) with qf and q'^ respectively. 
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Hence, the steepest-descent approximation naturally provides a mapping from the integration contour in the complex 
time plane onto a loop in the complex coordinate space connecting the endpoints [fig. [IJ. To avoid potential confusion 
with other methods, we emphasize that the appearance of complex paths has nothing to do with tunneling but rather 
is merely a consequence of the stationary phase approximation and holds also for systems with no barrier at all. In 
fact, it turns out that the complex semiclassical real-time trajectories used here never cross the barrier top, in contrast 
to paths emerging from ad hoc complexification procedures occasionally adopted to describe barrier penetration plj . 

Starting from the steepest descent conditions ( p"l| ) and exploiting Hamilton-Jacobi mechanics, one immediately 
derives 

Ps (0) = ip s (hl3) , p' s (0) = ip s (0) , E = E' = E. (12) 

Here, p s (u) \p' s (u)] is the momentum of the real time path q(u) [q'(u)] with energy E [E 1 ] connecting q s [q' s ] and qf [q'f]; 
accordingly, p s (o~) denotes the momentum of the imaginary time path q s (a) running from q' B to q s with Euclidian 
energy E = —p 2 /2M + V(q), see fig. [IJ Eq. ( p"2| ) can also be expressed as dT,/dt = with the solution 

^(Qf,l'fks,q' s ) = S(q f: q' f ). (13) 

Hence along the loop of steepest descent paths the full action is just given by the equilibrium action and thus 
independent of time. Differentiating Eq. ( [l3| ) with respect to qf,q'f, one finds 

p,(t)=ipo(h/3), p'M=iPo(0) (14) 

where po{o~) is now the momentum of the imaginary time path q~o(o~) connecting q'j with qf in imaginary time h(3. 
This path has Euclidean energy Ef which depends on qf,q'f and h(3 but not on t. Hence, we first deduce that the 
energies in Eq. (12) are given by Ef, which implies energy and momentum conservation throughout the loop in fig. [l]. 



Secondly, we arrive at the remarkable result that the sequence of time-dependent stationary phase points q s (t) [q' s (t)] 
is itself a minimal action path starting at q s (0) — qj [q' s (0) = q'^] with energy Ef. 

To complete the ordinary integrations in Eq. pi) over the initial coordinates qi,q[ we transform to fluctuations 
y = q% — q s and y' — q[ — q' s about the stationary phase points. An expansion of the full action ( [To| ) for fixed endpoints 
qf, q'f around the stationary phase points up to second order leads to £((?/, q'f\qt, q'j) = S(qf, q'f) + 5Y>( 2 >(y, y') with 

SWz(y,y')= 1 -(y,y')xW (15) 

where 

E(2) = / S SS E SS , \ 

is the matrix of second order derivatives, S ss = 9 2 S(<?i, q[)/dqf etc., to be taken at qi — q s , q[ = q' s . 

Inserting Eq. (^) into Eq. ([!]), the integrand now reduces to a product of Gaussian weight factors for deviations 
from the stationary phase points and an initial state factor 0(—q s — y)9(—q' s — y') describing deviations from thermal 
equilibrium at t = 0. Provided there is only one semiclassical path for each of the propagators, we obtain from Eq. (|l|) 
by virtue of Eqs. (|l^) and ([l5|) the semiclassical time dependent density matrix in the form 

pivfirffit) = ^P/3(qf,q'f)g(qf,q'f,t)- (17) 

Here, deviations from equilibrium are described by a "form factor" 

1 ru(qs) ru'(z : q' s ) 

g(q f ,q' f ,t) = - / dz / dz' e^ + z \ (18) 

where 



Det[gW] ; 



v/Dct[S(2)] 

with Det[E^ 2 ^] = E SS S S ' S ' — (E ss /) 2 . In deriving Eq. (Kjh we invoked that Hamilton Jacobi mechanics implies |2 



(19) 
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1/2 

= A{q f ,q' f ). (20) 

Note that for an initial equilibrium state, formally (?(•) — > 1 in Eq. (||) so that u,it' — > oo in Eq. (fl§l), the form factor 
becomes 1 and the semiclassical density matrix is in fact stationary. If there is more than one classical path one 
has to sum in Eq. (|l7]) over the contributions of all of them. Certainly, the above formulae ( |l8| ) and (|l9|) are only 
applicable as long as the Gaussian semiclassical and stationary phase approximations are valid, i.e. as fluctuations are 
sufficiently small. This will be seen to be no longer the case for low temperatures and/or very long times. How the 
classical paths in the complex plane can then be used as a skeleton for an extended semiclassical/stationary phase 
calculation will be shown below. 

In the remaining parts of the article we apply the general formalism to specific barrier potentials. Thereby, since we 
are particularly interested in the flux across the barrier, we restrict our investigation to non-diagonal end-coordinates 
qf and q'j — — <y close to the barrier top. This does not mean that we may constrain ourselves to study only local 
dynamics near the barrier top. Especially for lower temperatures, the nonequilibrium state in the barrier region is 
predominantly governed by global properties of the potential. 



A(q fl q s )A(q s ,q' s )A(q' fl q' s ) 
Det[£( 2 )] 



III. PARABOLIC BARRIER 



The semiclassical and the stationary phase approximations are always exact for quadratic potentials. Hence, as a 
simple test case we consider first a parabolic barrier 

V(q) = -\Mu? h q\ (21) 

Accordingly, the imaginary time dynamics runs in a harmonic oscillator potential. For the minimal action path q~o(a) 
connecting — qf with qf in time Hf3 one obtains 



This leads to the well-known equilibrium density matrix 

pp(qf , —q/) = — exp 
y47r<^ sin(cjbfif3) 



1 2 f 

- C ot(cj b ?il3/2)jL 



(23) 



with the relevant length scale 5b — y/h/2Mwb- 

The real-time dynamics simply follows. The classical real-time paths q(u) and q'(u) lead to the endpoints qf and 
—Qf, respectively, and hence obey q(t) = qf, q'(t) = —qf- On the other hand, the stationary phase condition ( |l2] ) 
implies q(t) = iq(Ti(3), q'(t) = iq(0) and we readily find q(u) = q(h/3 — it + iu), i.e., 

q ^ = ■ I \rio\ sin MW2 - it + iu)] , 

q'(u) =q(u + i%/3), 0<u<t. (24) 

At time t the imaginary time path qo(a) from —qf to qf is mapped onto the path q s (a) — q(ih(3 — ia),0 < a < K(3 
connecting q' s (t) = q'(0) with q s (t) = q(Q) (fig. ||). The stationary phase points 

q' s {t) = -q s (t)* (25) 

are as functions of t also classical paths moving away from the barrier top as t increases - q s (t) to the right and 
q' s (t) to the left for qf > 0. For longer times uJbt ^> 1, the stationary phase points asymptotically tend towards the 
limiting trajectories starting from qf — referred to as asymptotes henceforth. Similar to separatrices in classical 
phase space, these asymptotes divide the complex plane in regions of negative and positive Euclidian energy: in the 
sectors including the real axis classical real time motion has E < (but E < or E > 0) while in the remaining parts 
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E > 0. However, in contrast to simple classical separatrices the asymptotes are temperature dependent. The angel a 
of the q s and (^-asymptotes with the positive and negative real axis, respectively, is found as 



7T — LUbh/3 



(26) 



Now, for the nonequilibrium preparation ([|) the initial coordinates qi,q[ are constrained to Ke{qt}, Re{<^} < 0. 
Since q s (t) and q' s (t) are on different sides of the barrier, p(qf ,—qj,t) gains nonvanishing values only due to fluctuations 
that effectively shift qi away from q s and across the barrier top [see Eq. (|l9|)]. For the parabolic barrier potential the 
matrix elements in take the simple form 



Mw b [cot(uj b hp) - i coth(wfci)] , £ s v = S* s , 
Mujb 



so that the matrix X^ 2 ' can easily be diagonalized. One finds for the eigenvalues 



cot(ujbh(3) ± 



cot(uj b h(3) 2 : 



f 



sinh(o> hi)' 



1/2 



(27) 



(28) 



While, in principle, with Eq. ( |lq ) we can now evaluate the complete dynamics of the density matrix, we will concentrate 
here on the long time asymptotics of the nonequilibrium state. Then, in the asymptotic region uj^t 3> 1 the eigenvalue 
A_ tends to zero as A_ cx Mujb exp(— u>bt), reflecting the instability of the parabolic barrier. Hence, fluctuations around 
the stationary phase points with the least action increase occur in the direction of the eigenvector with eigenvalue 
A_. These fluctuations are of the form — |yi|exp[i(a + Ubhf3)] and y[ — \yi\exp(ia) so that qi and q[ move 
simultaneously along their asymptotes meeting at the barrier top (see fig. ||) . Now, inserting the matrix elements j27| ) 
and the stationary phase points ( |25| ) into Eq. ( JTgj ) and considering the limit u>bt 1, the relevant form factor turns 
out to be stationary 



\/7T 



i2qfl 



dx e 



(29) 



with n = v / cot(w 6 ?i/3/2)/(8 ( 5 j 
as 



The corresponding constant flux across the barrier is obtained from Eqs. (||) and (|T^ 



2ZM 



dq f 



q f =0 



which leads to the well-known result 



J B. = -r- 



4tt Zsm(L0 b h/3/2) 



(30) 



(31) 



Here, Z denotes an appropriate normalization constant which cannot be derived from the pure parabolic potential. 
This is not really a problem since realistic potentials always exhibit a well-behaved potential minimum and then 
Z follows e.g. as the relative normalization with respect to this minimum. Note that for the quadratic potential 
the results (|2^) and (|3l]) are formally valid for all times u>bt ^ 1- However, it was shown in (^] that due to the 
lack of a well-behaved ground state it makes physically only sense to use the parabolic barrier in T > T c where 
LOb^/k^T c — 7r. For lower temperatures Ubnf3 — > 7r large quantum fluctuations render the Gaussian approximation 
insufficient. Interestingly, the rate expression (|3l| ) diverges at the lower temperature T$ = T c /2 only where the 
parabolic density matrix (|23| ) ceases to exist. 



IV. DOUBLE WELL POTENTIAL 



A model well-behaved for the entire range of temperatures with many applications is the bistable dynamics of a 
particle moving in a double well potential 



G 



V(q) 



Mlo? , 



1 



(32) 



Here, the barrier is located at g = 0, the wells at q = ±q a , and the barrier height is Vf, = —V(q a ) = (Mu>f/4)q^. 
This potential exhibits rich quantum dynamics, namely, incoherent hopping between the wells over a broad range of 
temperatures that changes to coherent oscillations for T — ► 0. Due to the complexity of the dynamics this is a highly 
nontrivial problem for the semiclassical approach where the ratio Sb/q a serves as the small parameter. 



A. Thermal equilibrium 

The Euclidian mechanics in the inverted potential —V(q) can be solved exactly using Jacobian elliptic functions 
p5f . For the general solution one obtains 

q (q f ,a) = Bsn[oj(B)a- <p f \m], < a < hf3 (33) 

where the boundary conditions qo(qf,0) — —qf and q~o(qf.hf3) — qj fix the amplitude B and phase ipf. Since the 
potential is no longer purely quadratic - depending on temperature - there may be several solutions each of them 
with another amplitude. In Eq. (^3|) the frequency is given by 

B* 
29*' 



u(B)=uj b v/1-77 2 , V 2 = — , (34) 



and the phase can be represented as an incomplete elliptic integral 

r<if/ B i 

<t> f = F(q f /B\m) = / dx (35) 
Jo v (1 — str)(l — mr) 

with the so-called modul m = rj 2 /(l — rj 2 ). From the boundary condition q~o(hf3) = — qo(0) and the periodicity of 
the Jacobian function, sn[z + 2rK(m)\m] = (— l) r sn[z|m], r = 1,2,3,... with K(m) = F(l, m), the amplitude B is 
determined by 

uj{B)h(3 = 2rK{m) + [\ + {-l) r ]<t> s . (36) 

Since K(m), <pf > 0, for fixed uJb%(3 real solutions to this equation exist only for a finite number of integers r > 0. 

Let us briefly discuss the trajectories qo(qf, <j) as the temperature is lowered. For high temperatures only solutions of 
Eq. ( ^6|) with r = exist corresponding to direct paths from — qf to qf] particularly, (?o(0, cr) = 0. As the temperature 
drops below the critical temperature 

T c = hub/nkB, (37) 

i.e. LObhfi > ir, solutions of Eq. ( ^6|) with r = 1 arise. Then, for qf = the barrier top can be joined with itself also 
by two nonlocal paths denoted by q±(0,<r) oscillating in — V(q) to the right and to the left with amplitudes ±gi, 
respectively, and energy E\ = V(q\). With further decreasing temperature q\ grows and eventually saturates at q a 
for T — > 0. For finite qf the situation is rather similar: oscillating paths q±(qf,c) exist for all qf < q\. These paths 
connect — qf with qf via a turning point at ±31, thus, differing from q±(0,a) only by a phase shift. The described 
scenario repeats in an analog way at all T — T c /r, r = 2, 3, 4, . . ., where r counts the number of turning points. At zero 
temperature all these oscillating paths reach ±q a with the same energy E a = V(q a ) and are then called instantons. 

Now that all proper Euclidian trajectories are identified, the semiclassical equilibrium state follows readily. For high 
temperatures T > T c and end-coordinates qf near the barrier top, pp(qf,—qf) basically coincides with the parabolic 
result ( p3] ) and anharmonic corrections are negligible. This situation changes drastically for temperatures near T c . 
Then, the bifurcation of new classical paths leads to large quantum fluctuations and one has to go beyond the Gaussian 
approximation of the fluctuation integral. Slightly below T c a caustic appears for qf — q\. Since the region around 
T c was already studied in detail elsewhere p4| ] we omit here this crossover region and proceed with temperatures 
sufficiently below T c where near the barrier top Gaussian semiclassics is again applicable. It turns out that the paths 
newly emerging near T c are stable and dominate pp(q,q') for all T c > T > while the unstable "high temperature" 
path and those springing up at lower T give negligible contributions. Since for qf < q± all paths q±(qt, cr) differ only 
by a phase shift, one has for the corresponding actions 
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S ± (q f ,-q f ) = S+(0,0) = S-(0,0) 



(38) 



so that 



Pp(Qf,-Qf) = 2[A(q f ,-q f )] 1 / 2 exp[-S+(0,0)/h], 9/ < gi . 



(39) 



Accordingly, the matrix element ppfaf, — 5/) changes to a non-Gaussian distribution with a local minimum at qf = 
and two maxima at qf = ±qi. Thereby S+(0, 0) < 0, so that the probability ^(0,0) to find the particle at q = is 
substantially enhanced compared to its classical value. 

For T — > it is no longer sufficient to include only the trajectories with r = 1 in the semiclassical analysis but 
rather all other paths with r > 1 must also be taken into account. This is due to the fact that the smaller action 
factors of these latter paths are compensated for by zero mode phase factors from the corresponding fluctuation path 
integrals. Accordingly, all instanton contributions are summed up to yield e.g. for coordinates near the barrier top 



Pp(Qf,-Qf) 



-P[V(q a )+hu a /2]-Wa/h 



S a V2TT 



cosh 



cosh 



281 



sinh 



Here, W„ 



W(—q a , q a ) — —hf3E + S(—q a , q a ) is the short action for an instanton from — q a to q a . Further, 



cM-Wa/fr) 



(40) 



(41) 



denotes the WKB tunnel splitting with the well frequency iv a — ujb^/2 and 5^ = %/2Mw a . 



B. Dynamics of stationary phase points 

As in case of the parabolic barrier, the stationary real-time paths can be directly infered from the Euclidian dynamics 
at t = 0. From Eq. (|3^) and the stationary phase condition we have 

q,(t) = Bmtyf - iu{B)t\m], q' s (t) = -q s [(-l) r+1 t] (42) 

and qo(c) is mapped at time t onto q s {&) = qa[c + i(— l) r+1 t] where r follows from Eq. (|36|). In the sequel we 
always formulate the semiclassical theory in terms of the real time paths q s , q' s that "start" at the endpoints qf, —qf, 
respectively, and reach the initial points qi,q[ after time t. Since the endpoints <?/,— <7/ are fixed, while the most 
relevant initial coordinates depend on time, this backward view of the dynamics is in fact more transparent. The real 
time trajectories now start from the end coordinates we are interested in and lead to the relevant initial coordinates 
that need to be integrated over with an integrand weighted according to the initial deviations from equilibrium. The 
path q s (t) runs in the complex coordinate plane as a periodic orbit with period t p (qf) = 2K(1 — m)/uj(B) (fig. |^). 
Within one period it connects qf with qf via a loop crossing the real axis also after time t p (qf)/2 at the point 
q c (qf) — qs[qf ,t p (qf)/2] > q a . Thus q s (t) stays always on the same side of the barrier top and likewise q' s (t) on the 
other side so that the complex dynamics of the stationary phase points starting from qf and —qf, respectively, reflects 
a bounded motion in either of the potential wells. 

Let us consider the stationary orbits as the temperature decreases. For high temperatures T > T c , i.e. r — 0, each 
^/-dependent loop carries its own period t p (qf) and energy E(qf). If qf ^ 0, t p is small for T T c and the real time 
dynamics corresponds to a fast bouncing back and forth in the well. As the temperature is lowered the period grows 
while simultaneously the "width" of the loop q c {qf) shrinks. In the special case qf = the real time path reduces to 
a constant q s (fi, t) — 0. For temperatures T <T C the situation changes according to the appearence of new oscillating 
Euclidian paths q±(qf, <r) for qf < q\. In contrast to the high temperature case all stationary phase point paths with 
qf < q\ have then the same period t p (qt) = t p {q\) and energy E(qf) = E\, and differ only in their respective phases. 
Special cases are qf = and qj — q±: The path q s (qf, t), qf — > runs along the imaginary axis, while the orbit q s (q\,t) 
degenerates to a usual well oscillation along the real axis. These properties have a direct effect on the corresponding 
actions. One finds by employing Cauchy's theorem that after each period 

S[q s (qf,nt p ),q f ]= S[q s (qi,nt p ),qx], q f < q x . (43) 

Hence, all q s (qf,t) for qf < q\ can be seen as phase shifted copies of the specific real path q s (q\,t) having the same 
energy, period, and action increase during one period. In particular, the period t p {q{) is large for T<T C when q\ is 
still small, t p (q\) « ln(g a /gi)/wb, and drops down to t p (q a ) — 2ir /uj a in the limit T — > 0. 
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C. Nonequilibrium dynamics for high and moderately low temperatures 



For t — the density matrix is given by the initial state (||). The semiclassical time evolution of this state follows 
by inserting the proper classical paths into Eq.(|l7|). In the sequel, we mainly focus on the long time dynamics and 
are especially interested in a plateau region where the time evolution becomes quasi-stationary. 

We start by addressing the question when a plateau region does exist at all. Inserting the classical paths into 
g(qf,—qf,t) in Eq. (p^), the detailed analysis reveals that this function becomes stationary when the ratio q s (t)/p s (t) 
reduces to a constant. Since this will only occur within the parabolic barrier region, a least upper bound for a plateau 
region follows from the time interval within which a particle starting at a typical point qb near the barrier top continues 
to experience a nearly parabolic potential. This leads to t <C fp(<?&)- The lower bound is obvious: it is given by the 
transient time near the top, i.e. by l/u>b- Hence, a plateau region can be estimated to occur as long as a there is 
clear separation of time scales between local barrier motion and global well oscillations, i.e., 1 <C <-<Jbt <C u)bt p {qb)- 
Particularly, in the temperature domain where a plateau region exists, the Gaussian/stationary phase approximations 
are valid and we can actually calculate a rate. These approximations break down when the stationary phase points 
move far from the barrier top and times of order t p (qb) become relevant for the barrier crossing. This range will be 
addressed in the next section. 

For high temperatures T > T c the typical length scale qb can be identified with 8b = y / h/2Muib. Then, the 
separation of time scales fails for very high temperatures where k-QT>8Vb meaning that the thermal energy is of the 
same order or larger than the barrier height. With decreasing temperature t p grows so that for ojbhf3 of order 1 a wide 
plateau range appears with t p sa hi(q a /5b)/oJb- In the corresponding density matrix anharmonic corrections are small, 
and we obtain approximately the parabolic result (|29|). To get the rate, here, the proper normalization constant Z is 
taken as the partition function of the harmonic well oscillator 

Z = - e /3Vi> (44) 

2sinh(w o ft0/2) ' V ' 

Hence, from Eq. ( |30| ) one regains the known result 

= uJb_ sinhK^g/2) m 

2tt sm(Lu b h(3/2) { ' 

with the exponential Arrhenius factor and a characteristic ^-dependent prefactor that formally tends to uj a /ujb in the 
classical limit and describes the quantum enhancement of the rate as T c is approached. 

At this point we have to be very careful: a detailed analysis [ pTJ of the full density matrix p(qf,q'j,t), not only of 
its nondiagonal part, reveals that the nonequilibrium effects described by the flux state are restricted to the barrier 
region only in the presence of damping, consistent with the fact that finite temperature decay rates require coupling 
to a heat bath. In the absence of damping the full density matrix does not become quasistationary and the real time 
trajectories explore the strongly anharmonic range of the potential. Hence, an evaluation of the rate based upon a 
supposedly quasi-stationary flux state Pfl(<7/, — qf) for the undamped case corresponds to the transition state theory 
result. We refer to for a detailed discussion of this point. 

As the temperature reaches T c large quantum fluctuations occur and the impact of anharmonicities becomes sub- 
stantial. A detailed study of the bifurcation range around T c is quite tedious and was already presented in [ffijf . Thus, 
we omit explicit results here and proceed with temperatures T<T C where for coordinates close to the barrier top a 
Gaussian approximation - then around the newly emerging paths with amplitudes ±<ji - is again appropriate. As 
discussed above all real time paths with qf < q\ have now the same oscillation period t p {q±) f=a hi(q a /qi)/ujb. One 
observes that even though they are influenced by the anharmonicity of the potential via the Euclidian amplitude q\ , 
their time evolution for T<T C is still dominated by parabolic properties. Then, a somewhat lengthy algebra leads to 
the quasi-stationary density matrix 

Pftof,-<lf) = ^Pp(Qf,~Qf) 9 { n +) (qf,-<lf) + sjf ) 0/,-?/) (46) 

where describe the contributions from each of the two oscillating Euclidian paths. Note that due to symmetry in 
pp(qf, —qf), these contributions are identical and just lead to a factor of 2. In the temperature domain studied here, 
qi can be gained analytically from Eq. (pf) as 



Accordingly, one finds with some algebra for the thermal distribution 



P/»(0,0) 



exp 



mi 



i - 



(48) 



The form factor now has two contributions of the form ( |18|) and for K\ — d\o.{qi) / d{ij)bhjS) 3> 1 the corresponding 
integration boundaries read 



,(±: 



(?/) 



4(5 fcx /KT 



/ (±) (<j/, z) = m 



,(±! 



(?/) 



(49) 



Here, we employed that near T c the derivative E ss is dominated by d 2 S/dqf which is proportional to dE\/ d(u)bh(3) oc 
hl\. This way, using the normalization (44), the result for the rate is 



r = 



Mb sinh(w a ?i/3/2) 
2^ V2|sin(w6^yi 



-0V b 



(50) 



This expression is valid for temperatures T < T c where still «i ^ 1, a region which can be estimated as T somewhat 
larger than T c /2. There are two interesting observations to mention: first, the exponentially large term in the 
thermal distribution (|4^) - a consequence of the new Euclidian paths - is exactly canceled by a corresponding 
term which arises from the derivative of the form factor. This way, the rate is still dominated by the characteristic 
"Arrhenius factor" . Second, in the limit T — ► T c the above- T c formula ( f45| ) and the below- T c result (j50j ) both approach 
r c = (w&/27r) sinh(w 7i/3/2) exp(— /3V;,), however, the derivatives dT/dT are different. This discontinuity in the slope 
of the temperature dependent rate is removed by the full semiclassical theory (26) which takes the non-Gaussian 
fluctuations near T c into account and leads to a smooth changeover between the rate formulas ([45| ) and (^o|). 

With further decreasing temperature the amplitude q\ tends to saturate at q a so that n\ — > and the above 
rate expression is no longer applicable. Furthermore, the plateau region shrinks and eventually vanishes so that the 
assumption of a quasistationary flux state becomes inadequate even in the sense of a transition state theory limit 
of a weak damping theory. We note that in case of finite damping a meaningful rate can be found for much lower 
temperatures, then describing incoherent quantum tunneling. To investigate the time dependence of p(q/, —qf, t) with 
no damping in the limit of deep tunneling, we consider the case T = in the next section. 



D. Nonequilibrium dynamics for zero temperature 

Any Gaussian semiclassics to the real time propagator is expected to break down for very low temperatures and 
very long times where quantum tunneling comes into play. In fact, as yet no satisfactory semiclassical procedure was 
found to describe deep tunneling in the time domain. The crucial question thereby is: how can classical trajectories 
that either oscillate in one of the potential wells, here with energy E < 0, or move over the barrier, here with 
E > 0, produce exponentially small contributions to the semiclassical propagator which originate from quantum 
states connecting the two wells under the barrier. Here, we present a mechanism that is based only upon the complex 
plane mechanics discussed above and avoids any additional ad hoc insertion of "barrier paths" . Since the complex 
plane dynamics behaves as the usual classical real time mechanics, paths with E < never cross the barrier. However, 
a full semiclassical treatment needs to account for the dominant fluctuations about the semiclassical paths. Now, for 
T < T c there is a whole family of loop-like orbits in the complex plane; all with the same energy, period, and action 
increase after one period differing from each other only by their respective phases, i.e. by their crossing points qf < q\ 
with the real axis. It turns out that each time these orbits pass their end-coordinate qf there are other trajectories of 
this family arbitrarily close in phase space (see fig. ||) . The role of quantum mechanics then is to induce transitions 
between these orbits via small fluctuations. For sufficiently long times a path starting at a certain qf near the barrier 
top may successively slip down to an orbit with another phase q'f, eventually reach the stable regions around ±q a , 
and fluctuate in the long time limit between these regions. That this scenario actually describes the low temperature 
coherent tunneling dynamics has been discussed briefly in j2(| and will be described in some detail in the sequel. 

For T = the amplitude of the Euclidian time paths is q\ = q a . Thus, all stationary paths q s (qf,t), q' s (qf,t) have 
energy E — V(q a ) and period t a = t p (q a ) = lit juj a . Then, the Euclidian action S(qi, q^) suppresses energy fluctuations 
around E = V(q a ) exponentially, so that classical paths running in time t from qi ^ q s and q[ ^ q' s to qf and —qf, 

respectively, i.e. with E ^ V(q a ), are negligible. Further, studying the short action, W(q, q') = fjj dqp = S(q, q') + Et, 
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one finds according to Eq. ( (43[) that after each period W(qf,qf) — W(q a ,q a ) = 0. This result combined with the 
fluctuation prefactor [see Eq. (H)] gives for the Gaussian propagator after multiple round trips and for coordinates 
qf < q a 

\ G ^f^\ 2 « ntaiq L q} y " = 1.2,3,.... (51) 

Hence, the probability to return to the starting point decreases as the number of periods increases. In contrast, in 
the vicinity of the wells the Gaussian propagator coincides with the harmonic propagator. To be more precise, due to 
caustics in the semiclassics of this simple propagator at all nt a /2, an extended semiclassical analysis must be invoked 
leading to an Airy function; details of the procedure are well-known [ p7| and of no interest here. The important point 
is that in the barrier region the simple semiclassical return probability decays to zero for large times while in the well 
regions it remains constant. Thus, we conclude that the dominant quantum fluctuations neglected in the Gaussian 
approximation to the real time propagators are those that connect stationary paths with the same energy but different 
phases, i.e. initial coordinates qf. Effectively, these relevant fluctuations shift q slightly away from the classical path 
q s (qf,t) to reach another stationary path q s (q'pt) (cf. fig |3|). The corresponding change in action after a period and 
for small deviations is simply 

W(#,g/)«p a (g/,0)(#- 9/ ). (52) 

This repeats at subsequent oscillations. Hence, a "fluctuation path" can be characterized by its sequence of crossing 
points with the real axis after each round trip, e.g. by q^ ,k = 1, . . . , n for t — nt a where q^ — qj. Accordingly, 
a fluctuation path is not a classical path, i.e. it does not fulfill Newton's equation of motion, but can be seen as 
almost classical since it stays always in the close vicinity of a classical path. In the sequel we first explain the general 
structure of the extended semiclassical approximation and later on turn to details of the calculation. 

As an example let us consider a fluctuation path starting at qj that spirals around q a while the crossing point q^ 
with the real axis diffuses close to q a and returns to qf in t ^> t a (see fig. |]a). According to (|52| ) on the way to q a a 
particular path gathers an additional action W + (q a , qf) which is imaginary due to imaginary p s (q,0) [see Eq. dl3)] 
where 

\W+(q a , q f )\ - f 9 ' dq {2M[V(q) - V(q a )}} 1/2 (53) 

Jq f 

and the + [— ] sign stands for clockwise [anti-clockwise] rotation of the path in the complex plane. As long as the 
crossing point qf does not diffuse the imaginary time path connecting the endpoints of the two real time paths 
coincides at t = nt a with the imaginary time orbit connecting — qf with qf at t = 0. Taking into account the phase 
fluctuations, however, forces the endpoint of the imaginary time path to move with the endpoint q^ 1 of the "real 
time path" also towards q a - The mapped imaginary time path after t^>t a therefore runs from —qf to q a - According 



to (12) the additional amount of Euclidian action required for this deformation of the imaginary time path exactly 
counterbalances W(q a ,qf) so that the total action £ remains constant which reflects the stationarity of S along 
stationary paths. From close to q a the fluctuation path spirals back to qf. However, since q a is a branching point of 
the momentum there are two channels: the real time fluctuation path can maintain the direction of rotation or pass 
the turning point (TP) q a , thus changing the sense of rotation (cf. fig. [5^,). In the former case on the way back from 
q a to qf the fluctuation path crosses the real axis with the same direction of momentum as on the way to q a , so that 
due to W + (q, q r ) = — W + (q' , q) the path looses the action W + (q a , qf) again and returns to qf with W(qf,qf) = 0. In 
the latter case, momenta on the way back have opposite direction to those on the way forth so that the path arrives 
at qf with action W(qf,q/) — W + (q a ,qf) + W~(qf,q a ) = 2W + {q ai qf) and momentum — _p s (g/,0). 

Moreover, a fluctuation path starting at qf > can either move along the real axis to the right to reach q a 
or move to the left to arrive at — q a . In the latter case, the crossing point q^ diffuses across the barrier top so 
that the path initially spiraling around q a finally orbits around — q a with opposite sense of rotation. Accordingly, 
due to W + (q a ,0) — — W~{— q a ,0) the real time action factor exp[iW(±q a , 0)/%] grows or decreases exponentially 
for diffusion to the right or to the left, respectively. In any case, near ±q a the semiclassical propagator has to 
match onto the propagators in the harmonic wells. For the two lowest lying eigenstates which are relevant here, the 
matching procedure was discussed in detail by Coleman |p8[ . Correspondingly, these two states determine the relevant 
propagator in its spectral representation. Then, it turns out that a TP may only occur if iW(qf, ±<? a ) < 0. This has 
profound consequences on the extended semiclassical approximation: (i) A relevant fluctuation path from qf to qf 
must reach ±q a rotating clockwise to have a TP. Then W + (q a , qf)+W~(qf, q a ) — 2i\W(q a , <7/)| and the corresponding 
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contribution to the propagator has an exponentially small factor exp[— 2|VF(<7 a , 5/)|]- (ii) A fluctuation path with more 
than one TP has to alternatively visit TPs at ±g a thereby changing its sense of rotation repeatedly. In passing from one 
TP to the next the path gathers the action W (0, -q a ) + W + (q a , 0) = W~(0, q a ) + W + (-q a , 0) = i\W(q a ,-q a )\ = iW a 
which coincides with the instanton action in A introduced in Eq. ([4l|). Hence, according to (i) contributions from 
fluctuation paths with TPs do not play any role for short times. For longer times, however, they become increasingly 
important, particularly, since a fluctuation path may spend an arbitrary period of time at the TPs ±q a where 
V'(q a ) — before leaving them. The detailed analysis shows (see below) that for t 3> t a the phase space of equivalent 
fluctuations with one TP is therefore oc t which compensates for the exponentially small action factor. Moreover, at 
each TP a path gathers an additional Maslov index v — > Then, according to (ii) the full density matrix is given 

by a sum over v, v taking into account the proper order of TPs, i.e. 

P(<lf,-<lf,t) = Pw'faf'-Qf't)' ( 54 ) 

i/,i/'>0 

where p v y denotes the contribution from relevant fluctuation paths with v TPs in the forward and v 1 TPs in the 
backward propagator. Note that v 1 v' only label the TPs of the two real time paths. For each v, v' one has also to 
sum over all imaginary time paths connecting the endpoints of the real time orbits. 

To evaluate the sum in Eq. ( |54| ) we start by analyzing the term with v = v' = 0. As discussed above the diffusion 
of the real time orbits is then irrelevant and the imaginary time path has to run form — qf to qf. To lowest order 
in A [see Eq. p0|)] we have the two imaginary time paths q±(qf,c) emerging at T c that connect —qf with qf via 
TPs at ±q a , respectively. At T = other solutions of the imaginary time dynamics with r > 1 [cf. Eq. (|36|)] just 
contain additional intermediate instantons, i.e. imaginary time trajectories connecting q a with —q a or vice versa. 
For an equilibrium initial preparation (small) fluctuations about the stationary paths towards ±q a , respectively, give 
identical contributions and we thus recover pp(qf, —qf)/Z [Eq. (^|)] in the limit T — > with the partition function 

Z = 2exp[-(3V(q a ) - Phu a /2] cosh(fy3A/2). (55) 

For the nonequilibrium preparation, however, only fluctuations towards — q a contribute (fig. |^b). This way, one obtains 
for coordinates qf near the barrier top 



^Mlf) 2 (56) 

where we used the short action of an instanton W a = W(q ai —q a ) with W a = \W(q a ,— q a )\ and ipo(qf) denotes the 
semiclassical ground state wave function in the double well. 

The next order real time paths are those with one TP, i.e. v = 1, v' = and v = 0, v' = 1 in Eq. (|54|). Thereby, 
the real time path q s (t) makes an excursion from qf via a TP at q a to — q a in case where at t = the endpoints — qf 
and qf are connected by an imaginary time path <f_, while it diffuses from qf via a TP at —q a to q a in case of 
Accordingly, one observes that for an equilibrium initial preparation all contributions cancel, e.g. the contribution 
corresponding to g_ with v = 1, v' = cancels that corresponding to q+ with v = 0, z/ = 1. In fact, it can be 
shown in the same way that for an equilibrium initial state all terms in the sum (|54|) with u, v 1 > vanish. However, 
for the initial preparation (j^) a finite result follows due to the projection onto the left side of the complex plane. 
Hence, both real time orbits have to end near — q a whereby one trajectory has a TP at q a . According to the above 
discussion we gain the following action factors: For v = 1, v' = one has exp[— 3|1 / F((7 a , 0)\/h + q a qf / (46%)] from the 
forward and exp[~\W(q a ,0)\/h + q a qf / (AS^)] from the backward propagator (cf. fig. ||c) while v = 0, v' = 1 gives 
exp[-\W(q a) 0)\/n - q a q f /(ASl)} and exp[-3\W(q a , 0)\/h - q a q } /(A& 2 a )] (cf. fig. |d), respectively. After expanding the 
integrand in (|lj) around — q a up to second order, the ordinary integrations over the initial coordiates are seen to be 
restricted to the harmonic range around — q a so that the 8 functions can be put to 1. Then, the integrals effectively 
describe the stationary real time motion of the equilibrium well distribution pp(—q a , —q a ). Combining these findings 
yields 

pi{<H>-Qf,t) = pi, (qf,-qf,t) + p ,i(g/, -?/,*) 

= 18 *(t) ex V [-2\W(q a , -q a )\/h] smh(q fqa /2S 2 a ) p f3 (-q a , -q a )/Z. (57) 

Here pp(—q a , ~q a ) includes a sum over multi-instanton contributions of the imaginary time paths, resulting for T — > 
in pp(—q a , —qa)/Z — y/ 'Mu) a / 'AttTi. Further, $(<) oc t takes into account the phase space contribution from equivalent 



Po,o(qf,-Qf,t) 



j™ yPp(lf,-Qf) 

Z p—>oa Zj 



2nS a 



exp[-W a /h] cosh (q f q a /A6l) 
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fluctuation paths connecting qj via a TP at q a with — q a . These paths differ only in their "sojourn times" at the TP 
q a . To evaluate $ we adopt the method outlined in p^] to which we also refer for further details and write 

G t (q a) -q a ) = I duG t ^ u {q a ,Q)q s (Q,u)G u (0,-q a ). (58) 
Jo 

Since for t 3> t a the sojourn time is exponentially large, one can actually sum the intermediate time step u over the 
entire time interval up to negligible corrections. To calculate the semiclassical propagators in the integrand of d58| ) 
one exploits that a fluctuation path moving from q = to q a in time t 3> t a spends almost all time by orbiting in 
the vicinity of q a thereby diffusing along the classical stationary paths ( ff2| ) with q/^Sq a towards the TP. Hence, the 
time dependence of the propagators is determined only by the asymptotic behavior of the fluctuation paths. Then, as 
already derived in the previous paragraph, the actions turn out to be independent of time up to exponentially small 
corrections and their sum gives rise to the action factor exp[— ^(gq, — q a )\/h] in Gt( g a , —q a )- For the prefactors one 
uses the representation (^), and exploits the fact that g s (0, u)^ A t - U (q a , 0)A U (0, —q a ) depends for t 3> t a on time and 
temperature only through exp[— (ft/3 + it)w a /2] while its dependence on the intermediate time step u is exponentially 
small. This way, since both exponential factors are already accounted for in (|57|), we arrive at 

$(t) = G t (q a , -q a ) exp[\W(q a , -q a )\/h + (ft/3 + it)w«/2] (59) 

which leads to 

*(t)=t*£*L. (60) 



Combining this result with Eq. ( |57| ) we derive the one-TP contribution to Eq. (]54|) as 

Pl(q f ,-q f , t) = itA 2 smh{q a q f /2S 2 a )/(uj a q a ) (61) 

where the tunnel splitting is specified in (|4l|). Likewise, contributions from real time path with more than one TP 
can be calculated where the proper order of TPs must be taken into account. Eventually, only contributions with 
i 2k ~ 1 A 2k , k = 1, 2, . . . survive and the time dependent density matrix ( [54] ) for coordinates near barrier top reads 

p(<lf,-qf,t) = -ipo(q f ,-q f ) 2 + iA sin(At) smh(q a q f /25 2 a ) /{u a q a ). (62) 

Hence, the initial state (0) develops an imaginary, time dependent part from which the tunneling current (J^j is gained 
as 

J(t) = Asin(At) (63) 

describing coherent tunneling between the wells. This shows that a systematic semiclassical analysis of the real time 
dynamics of the system covers also low temperature tunneling. We note that the tunnel splitting coincides exactly 



with the result of the instanton approach [see Eq. ( 41 )] where A is related to the action of an imaginary time path. 
Within the real time description the "instanton dynamics" is replaced by the above-mentioned diffusion along the 
real axis in the complex coordinate plane. 

Before we conclude this section let us briefly sketch how the crossover from coherent decay to incoherent tunneling 
occurs within the present formalism as the temperature is raised; details of the calculation will be presented elsewhere. 
For finite temperatures T > the energy of the stationary paths is \ — |^(<7i)| < |V(g a )| so that the TPs of the 
fluctuation paths q\ < q a are shifted towards the barrier top. Accordingly, V'{q\) ^ and the corresponding actions 
in Eq. (|58| ) are no longer independent of time. Thus, the time interval t S> t a may eventually exceed the region where 
the integrand gives a contribution so that <I>(i) — > $a saturates. We note that for a system with dissipation basically 
the same mechanism, namely, an effective action depending on time via damping induced correlations, may cause 
incoherent tunneling even at T = 0. 

V. ECKART BARRIER 

As another instructive example we analyze in the following the transport across a genuine scattering potential, 
namely, the so-called Eckart barrier 
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Here, Vq is the barrier height and Lq the typical interaction range. We drop the condition V{q = 0) = in this section 
so that energies are shifted by Vq. In fact, the real time dynamics in this potential is much simpler as in the double 
well: particles steadily injected from a thermal reservoir to the left of the barrier built up a flux across the barrier 
that is stationary for all times after a certain transient time has elapsed. Thus, the corresponding quantum dynamics 
is described by a barrier transmission rate for all temperatures. In a semiclassical expansion we use h/ \/2MLqVo as 
the small parameter which demands high and broad barriers. 

A. Thermal equilibrium and stationary phase points 

The solution of Newton's equation of motion for the Eckart barrier in imaginary time reads p0[ 



q (q f ,a) = L arsinh | \j V ° g — sin [uj{E)u - 4> f ] j (65) 
where we introduced the energy dependent frequency 

/ K 

(66) 




with the barrier frequency u>b = -\/2Vo/MLq. Energy E and phase <f>f are determined by the boundary conditions 
go(<7/)0) = — <7/ and q~o(qf,hf3) — <//. Accordingly, employing one of these conditions to fix <f>f, the energy can be 
evaluated from 

w(E)hp = rir + [1 + (-l) r ]^/, (67) 

where for given temperature real solutions exist only for a finite number of integers r > 0. Accordingly, this nonlinear 
equation gives the amplitude of the semiclassical path. A s a fu nction of temperature the solutions ( |65| ) and (|67j ) 



resemble those in the inverted double well potential (cf. Sec. [IV A| ). Particularly, for T < T c again all paths with r = 1 
have the same energy independent of qj 

Ei ^ V( qi ) = J^* (68) 

with amplitudes ±q± and the same frequency io\ = u){E\). An important difference to the double well case, however, 
is that the amplitude q\ of the paths here grows without any limit as T — > 0, i.e. E\ — > 0. Hence, the equilibrium 
density matrices pp{qf , —Qf) differ qualitatively in the deep tunneling regime. While for the double well potential 
near T = contributions from all multi-instanton paths must be summed up [see Eq. (|40|)], here, the density matrix 
is dominated by the oscillating paths newly emerging around T c = uj^h/k-Q-K for all T < T c . For further details of the 
Euclidian semiclassics we refer to J3(i| . 

Now, for the stationary phase points we obtain 



q s (t) = L arsinh { J -^r~ sin [<Pf - iu(E)t] , q'M = ~q s [(~l) r+1 t] (69) 




which are connected by the Euclidian path q s (<r) = qo[(J + i(— l) r+1 i\. Starting at q/ [—q/] the path q s (t) [q' s (t)] 
describes for large times an almost free motion parallel to the real axis in accordance with the asymptotically vanishing 
interaction V(q) — > as q — > ±oo. The energy of q s (t) is controlled by temperature where qualitatively the two ranges 
T > T c and T < T c must be distinguished. In the first case, E depends on qf and for gj > we find asymptotically, 
i.e. for u>(E)t > 1, 



q s (qf,t) « L arsinh 



sinh( g/ /£ ) (@)t 
2sm{uj b fiP/2) 



'tt uj b hf3\ 

-^ [- — ) (70) 
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so that for T > T c all stationary paths are restricted to the strip i(Lo/2)[— (tt — u>bh(3), (it — u)bh(3)] (cf. fig. ^). In the 
range T < T c and for end-coordinates qf < qi the oscillating Euclidian paths determine the energy independent of 
qf as E = E\ [see Eq. jS8|)] with the frequency u>i = u>(Ei). Hence, in T < T c all path starting in the barrier range 
{qj<Lo) nearly coincide asymptotically (see also fig. ^) for uj^t ^> 1 

q s (q/,t) faqi+L In [sinh(wit)] + iL n/2. (71) 

Orbits with qf > Lq come also close to ILq-k jl but on a larger time scale, e.g. for qf close to qi 3> Lq on the scale 
1/lui. Since Im{^} — > —iLQit/2 for large t we conclude that the real time stationary phase dynamics for T < T c is 
restricted to the strip [— iLq-k/2, iLq-k/2] in the complex plane (fig. ||). To complete this discussion we address the 
special cases qf = qi and qf = 0, respectively. In the former case the motion starts with zero momentum and takes 
place along the real axis, asymptotically (u>it ^S> 1) as 

q s (t, qi) « qi + L u>it. (72) 

In the latter case the orbit can only be defined as the limiting trajectory of q s (t, qf) for qf — > 0, thus, running along 
the imaginary axis from q = to q = iLq-k /2 and afterwards parallel to the real axis. 

The described complex plane dynamics depends essentially on the analytic properties of the potential V(q). Inter- 
estingly, in case of the Eckart barrier one has for complex q the periodicity 

V(q) = V(q + iL nir), n integer. (73) 

Hence, the complex plane falls into strips [(2n — l)iLoir/2, (2n + l)iLoir/2],n integer, parallel to the real axis each 
of which with identical classical mechanics and corresponding stationary phase paths q s (qf + iLomr,t). As shown 
above, for T < T c the real time dynamics starting from the real axis at t = reaches asymptotically the boundaries 
of the strip n = 0, while in the range T > T c it does not. This has crucial impact on the semiclassical analysis for low 



temperatures as will be discussed in Sec. V C 



B. Stationary flux for high and moderate low temperatures 



The systems starts from an initial state where the thermal equilibrium is restricted to the left of the barrier, thus 
extending to q — > — oo. Accordingly, after a certain transient time has elapsed the flux across the barrier remains 
stationary forever. Since with increasing time the stationary phase points move away from the barrier top, for large 
times fluctuations of the order of Lq or larger are needed to shift qi into the region Re{qi} < 0, then rendering the 
Gaussian stationary phase approximation insufficient. Yet, we can use Eq. ( |l7| ) to gain Jr as long as the flux becomes 
stationary on a time scale within which q s remains smaller than Lq. 

In the range T > T c and for qf near the barrier top one has E rj Vq so that u>(E) ss a>&. Then, the density matrix 
p(qf, —qf, t) tends to stationarity on the scale 1/wf, while q s (t) reaches Lq on the much longer time scale ]n\Yo/(E— Vq)] 
only. We thus regain within this time window approximatively the parabolic result ( p9|) in the semiclassical limit 
large Vq and Lq - where anharmonicities are negligibly small. Correspondingly, the rate reads as specified in Eq. ( |3~l| ) 
with ui b = V / 2V /MLl. 

For temperatures T < T c the transient time range grows according to l/tui = hf3/ir, while the upper bound for the 
validity of the Gaussian approximation eventually shrinks to 1/uib- Hence, while one can no longer use the local barrier 
dynamics for temperatures T T c , in a region sufficiently close to T c a rate calculation along the lines described in 
Sec. IVC still makes sense. Accordingly, for T<T C the density matrix is obtained as in Eq. ( ff6| ) with the amplitude 
qi derived from Eq. (B7j) for r = 1 as 



qi = Lq arsinh 



This way, one gets the rate 



(74) 



Jfb_ 

\-K 2 Z 



arsinh y {ubh(3 / tt) 2 



1 



exp(-/3F ) 



(75) 



for temperatures below T c but still above T c /2. For even lower temperatures higher order terms in the expansion 
around the stationary phase points q s ,q' s must be taken into account. In the following section we show that for T 
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below T c /2 the rate is dominated by quantum tunneling which requires an extended semiclassical analysis. Thus, 
a higher order expansion is needed only in the close vicinity of T c /2 where the changeover from the thermal to the 
quantum rate occurs. 



C. Stationary flux for low temperatures 



The breakdown of the Gaussian stationary phase approximation for lower temperatures indicates also a breakdown 
of the simple semiclassical approximation to the real time propagators for large times. In fact, one needs to carefully 
analyze the quantum fluctuations around the classical paths to capture tunneling processes. In the sequel we proceed 



in the spirit of Sec. [VP and search for relevant phase fluctuations. 

We begin by recalling that in the range T<T C and for coordinates qf < q\ the classical mechanics in the complex 
plane takes place in strips [(2n— \)iLq-k j2, (2n+\)iL^n /2], n integer, parallel to the real axis. One thus has families of 
classical paths (cf. fig. ||) all with the same energy E\ that start at t — > — oo to the far right on the lines (2n — l)iLoir/2, 
run close together with almost vanishing momentum — MLqloi towards the barrier top, pass at t = the coordinates 
qf + inLgTr and then leave again to the far right moving close together with momentum MLqU>i asymptotically along 
the lines (2n + l)iLo 7r /2. Accordingly, for T — > in classical phase space, see fig. 0, orbits with different phases 0/, 
i.e. different qf, but from the same or from adjacent strips lie arbitrarily close to each other in the asymptotic range 
where |I^(<z)| — > 0. The effort of quantum fluctuations then is to link these paths which reflects the asymptotically free 



particle diffusion in the Eckart potential. In simple semiclassical approximation one has asymptotically the propagator 

/M/^-^ (76) 



M Y" r.„(g-gO 



Gt(9 ' 9 ' )= <2^) ° xp 



2H 



so that for fixed q — q' and large times the transition probability decreases as \G t (q,q')\ 2 cx 1/t [cf. Eq. (jsil)]. 
Correspondingly, two different types of fluctuations can be identified: one type of fluctuations connects paths q s (qf + 
iLoniT,t) and q s {q'f + iLQnn,t) within the same strip, while the other type of fluctuations switches between paths 
1s(qf + iLonir,t) and q s (q'f + iL<j(n + l)n,t) in adjacent strips. The first type is already accounted for in the simple 
semiclassical approximation to the real time propagator since these fluctuations never leave strip n = and stay in 
the close vicinity to the asymptotic q s (qf, t). In contrast, the second type is relevant beyond the Gaussian semiclassics 
since it causes large deviations and allows a path q s (qf, t) by subsequently diffusing to another strip to reach a path 
q s {q'f + iLoJiiT,t) with q'f far from qf and n large. Interestingly, this second kind of fluctuations does not exist for 
T > T c where asymptotically there is always a gap iuiffifi between paths in adjacent strips (see fig. ^jj). 

As an example let us consider a trajectory q s (qf,t) with qf close to the barrier top for T <C T c and times t 3> 
l/u)\ 3> l/wft. In this limit the orbit runs for wyt 3> 1 along the boundary iLq-k /2 of the strip n — where fluctuations 
of the second class bridge the tiny gap to an orbit q s (q'f +iLoTr,t) with a different q'f in the strip n = 1. This trajectory 
passes q'f + iLqtt, and exploiting the periodicity of V(q) the corresponding change in action W(qf,q'f) is shown to 
read as in Eq. (^). Obviously, the described fluctuations always lead from an outgoing to an ingoing orbit thereby 
increasing the strip number which in turn requires a momentum fluctuation of order 2|<j s (t)| = 2MLquJ\. Estimating 
typical momentum fluctuations by Tl/Lq one rederives from %/Lq ^> MLqlui the condition T <C T c so that at low 
temperatures these fluctuations will indeed occur. By the same procedure the path q s {q'f +iLoir, t) can be linked to a 
path q s (q'f +iLo2Tr,t) and so forth. Similar as in case of the double well potential a "fluctuation path" is characterized 
by its sequence of crossing points q^ + iL kn, k = 0, 1, 2, . . . , n [q^ — qf], with the lines iL kn, i.e. the copies of 
the real axis in the strips k. Accordingly, for very large times t 1/uj\ the point -\- ILqut: moves with increasing 
n along the positive imaginary axis while simultaneously q^ can slide down the real axis and away from the barrier 
top to reach the proximity of q\. From close to qi relevant fluctuation paths traverse q\ as a turning point (TP) - q\ 
is a branch point for the momentum - and return via the described scenario to qf in the strip k = 0, however, crossing 
the lines iLokir with opposite direction of momentum as on the way forth. The total change in action is imaginary 
and given by W(qf, qf) — W(q\,qf) — W(qf,q\) = 2W(q\,qf), n arbitrary but large, where for qf close to the top 

\W(qi,q f )\ = dq[2M(V(q)-V( qi ))} 1/2 

u b ML qf. (77) 



Jqj 

ttVq 

UJb V U>bfof3 



In a similar way, the sequence of q^ of a fluctuation path starting at qf can move directly towards the barrier top, 
diffuse across the barrier to enter the left halfplanc of the complex plane, and end up in the asymptotic region Rejq} — > 



16 



— oo. Since in leading order the semiclassical propagator has asymptotically to match onto the free propagator fl76|), 
a TP may only occur if iW(qf,±q±) < 0. Hence, what we discussed in Sec. IV D [see paragraph above Eq. (M)] can 
directly be transfered to the situation here and the density matrix can be cast into the same form as in Eq. (|54|). In a 
notable difference to the double well potential, however, the TP qi here is not an isolated extremum of the potential 
meaning that each TP - for Euclidian and real time fluctuation paths as well - is not related to an additional phase 
factor for equivalent paths. 

After having elucidated the general structure of the semiclassical density matrix we now turn to the explicit 
calculation of the sum (54) and begin with the term po,o(<Z/j — If, t)- This matrix element follows by the same arguments 
as given in Sec. IV D. Since the equilibrium density matrix for the Eckart barrier is dominated by the oscillating 



paths newly emerging around T c for all lower temperatures, all further contributions from Euclidian trajectories with 
r 7^ 1 in Eq. (p7\) are negligible. Accordingly, we find for coordinates around the barrier top 



Po,Q{Qf>-Qf) 



i f 



ttV 



1/2 



zl [ oj 2 b n 2 p[(co b n(3y(i ~ q }/L%) - ^ 



exp 



ttVq 



(Vbn/3 



(78) 



Note that in contrast to bounded systems the above density matrix remains temperature dependent even for T <C T c . 

To next order real time paths with v — l,v' = and v = 0, v' = 1, respectively, contribute (cf. figs. ||c,d). For 
t l/u>\ a relevant real time fluctuation path with v = 1 starting at g/ moves via a TP at q\ + iLon'ir, ml large, to 
the left halfplane where it eventually crosses — q% +iLonn, n large, to run along the line iLoiiir and reach qi +iLomr in 
the far left. For the segment of the fluctuation path from qj via a TP to —qi +iLonw the corresponding action factor 
is exp[—3\W(qi,0)\/h + uji,MLoqf/h}. Due to the periodicity ( f73| ) of the potential the segment from —q\ + iLqu-k to 
qi + iLgnn, qi < —qi, can just be treated as the corresponding one along the real axis; for very large times io\t 3> 1, 
i.e. \qi\ 3> qi according to Eq. (f72|), we then get the action factor exp[— iMq 2 /2ht}. Hence, the corresponding relevant 
real time propagator reads 



G t (qf,qt) 



A{q f ,qi) exp 



3\W(q u 0)\ 



"bMLpgf q[_ 

h 2ht 



(79) 



Similar, the propagator from —qj directly to q[ + ih^n-K is gained. The crucial point is now that for the integral in 
Eq. (Q) there are no longer isolated stationary phase points but rather all qi , q\ on the line iL^nir and to the far left of 
the barrier top make the integrand for very large times stationary. The ordinary integrals in Eq. ([j]) can thus be seen 
as sums over stationary phase points qi , q[ whereby their distance is weighted by the asymptotic thermal distribution, 
i.e. in leading order the free particle equilibrium density matrix 



pp(q,q') = 



M 



2irh 2 p 



1/2 



exp 



M{q-q'f 



2h 2 (5 



(80) 



Accordingly, for T — > one has Pf3(qi,qi) — > S(qi — q'i) so that contributions from qi ^ q\ are caused by thermal 
fluctuations at elevated temperatures. Further, for large qi, q[ and large times the prefactors A(qj , qi) and A'(—qf, q^), 
respectively, are independent of qi,q' iy thus allowing us to carry out the qi,q'i integrals over the exponentials only. 
Then, using —q\ as an upper bound for the asymptotic coordinate range it turns out that for u>xt 1 the result for 
the integrals in leading order is irht/M. Now, combining all factors we finally obtain the time independent density 



Pl(Qf,-Qf) 



lim p 1) o(qf,-q f ,t) + po,i(qf,-qf,t) 

uif>l 



ZLn 



47rVb 



1/2 



smh(2ujbMLoqf/K) e 



-4|W(gi,0)|/fi 



(81) 



where |W(gi,0)| follows from Eq. ([77]). Employing the same procedure, contributions in the sum ( pi} ) from real 
time paths with more than one TP can be derived, however, they contain additional action factors and are thus 
exponentially small compared to p%. Hence, the stationary semiclassical density matrix for low temperatures and very 
large times is found as 



Pa(qf,-Qf) 



2Z 



pp(<if,-qf) + pi(qf,-qf) 



(82) 
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with p/3 as specified in Eq. (|r8|). Finally, from Eq. (^) we gain the thermal tunneling rate 



1/2 



-i\w{ qi ,o)\/h 



(83) 



This simple formula is applicable as long as q\ > Lq, a temperature range which can be estimated by T below 
T c /2, or equivalently ojbftfi > 2n. To be precise, there is also a lower bound for the temperature. Namely, for 
T — ► any semiclassics in the Eckart barrier breaks down due to the fact that then tunneling takes place in the low 
energy range near the base of the barrier where the wave length of a wave function tends to exceed the width of the 
barrier. From the known exact transition probability [see e.g. ]30| ] one derives that this scenario becomes relevant 
for aji,h/3 3> 2n (Vo/hu>b) corresponding in the semiclassical limit Vo/fttUb 3> 1 to extremely low temperatures. In the 
broad temperature range between these bounds, i.e. 2tt < uibfi-ft^^ir (Vo/hoJb), the above rate expression describes 
the decay rate with remarkable accuracy when compared to the exact result even for moderate barrier heights (see 
fig. ^). Table I presents a numerical comparison with results from other approaches. For temperatures above T c /2 the 
real-time semiclassical rate is slightly too small and coincides for T > T c with the well-known "unified" semiclassical 
rate formula gained by the thermal average over the transmission T(E) — 1/{1 + exp[S '(E) / ft]} , where S(E) is the 
bounce action for T < T c /2. The small deviations from the exact rate are due to the fact that in the simple version 
of the theory presented here anharmonicitics of the potential are neglected for T > T c and taken into account only 
in leading order in T c > T > T c /2. A perturbative expansion in an anharmonicity parameter allows for a systematic 
improvement. For the same reason, the temperature region around T c /2 is not well described. In the deep tunneling 
region T < T c /2, which is notoriously problematic for real-time rate theories, our theory performs excellently. In fact, 
the low temperature formula ( |S3| ) turns out to be identical to the result derived within the instanton/bounce approach 
sec e.g. @]. The bounce is an oscillating Euclidian orbit, periodic in phase space, which connects qf — with itself, 
thus emerging as a solution of Eq. @ at T = T c /2 (r = 2). While in imaginary time methods the bounce trajectory 
describes barrier penetration, here, effectively the same tunneling rate arises from fluctuations around real time paths 
the energy of which is fixed by oscillating Euclidian orbits, closed in phase space, emerging at T c (r = 1). These 
latter minimal action paths solely determine the semiclassical thermal equilibrium for lower temperatures T < T c , 
thus establishing within a semiclassical real time approach the relation between the thermal density matrix and the 
thermal tunneling rate, since long an open question in thermal rate theory. 



TABLE I. Transmission factor P 
same as in fig. ]A 



r/r c ; for the symmetric Eckart barrier. T c i is the classical rate and parameters are the 



UJbftP 


pa 


P h ■ 

' uni 


■Pqtst 


"SQTST 


pc 
± ex 


1.5 


1.10 


1.10 


1.13 


1.13 


1.13 


3 


1.50 


1.50 


1.54 


1.52 


1.52 


5 


2.98 


3.84 


3.18 




3.11 


6 


3.86 


8.92 


5.74 


2.2 


5.2 


8 


21.99 


17.97 


29.3 


11.9 


21.8 


10 


136.2 


132.2 


248 


149 


162 


12 


1613 


1606 


3058 


3006 


1970 


16 


6.03 • 10 s 


6.03 • 10 5 




2.56 • 10 6 


7.41 • 10 5 


18 


1.54 • 10 7 


1.54 ■ 10 7 




9.1 • 10 7 


1.88- 10 7 


a p 
1 rsemi 


is the transmission factor 


as derived by the real-time semiclassical approach presented 


in this paper 





b ^uni is the transmission factor of the "unified" semiclassical approach 

c Pqtst is the transmission factor according to the simplest version of Pollak's QTST, from Ref. 

^Psqtst is the transmission factor according to the full semiclassical version of Pollak's QTST, from Rcf. j32j 

e P C x is the exact transmission factor 



18 



This represents substantial progress when compared with other attempts. While Pollak's new quantum transition 
state theory [Tl| is based on a numerically exact evaluation of the thermal flux, it suffers from a simple semiclassical 
approximation to the real time propagators. In the associated full semiclassical calculation by Pollak and Eckhardt p3] 
only half of the bounce action appears in the exponential factor for temperatures below T c /2 and the corresponding 
tunneling rates are too large (see also Table I). The centroid method gives the correct action factor, however, 
its semi-empirical factorization of thermal and dynamical contributions leads to a prefactor which is too small for 
lower temperatures. Finally, from semiclassical real time calculations for the Eckart barrier based on the simple 
semiclassical propagator ]T^ , p^| tunneling probabilities in the deep tunneling regime cannot be properly extracted 
since they depend strongly on the initial state. 



VI. CONCLUSIONS 



We have developed a unified semiclassical theory that describes the real time dynamics of quantum statistical 
systems for all temperatures including coherent and incoherent processes. Starting from the exact nonequilibrium 
dynamics the approximate density matrix is gained by employing semiclassical propagators in real and imaginary time 
combined with a stationary phase evaluation. Accordingly, the relevant classical mechanics takes place in the complex 
coordinate plane where the energy of classical real time trajectories is fixed by the Euclidian orbits determining the 
equilibrium distribution. Consequently, real time paths follow by solving an initial value rather than a boundary value 
problem. While this procedure can be used to study the dynamics for a wide class of systems and initial conditions, 
here, we concentrated on the flux across a double well potential and an Eckart barrier. Then, for high to moderate 
temperatures the Gaussian approximation suffices to obtain a stationary flux. In the tunneling domain, however, this 
approximation fails and the complex plane dynamics allows to identify the dominant quantum fluctuations in the 
real time propagators. These are zero-mode like phase fluctuations which give rise to a diffusion along the scaffold 
of classical orbits. Quantum tunneling in the real time domain can semiclassically thus be interpreted as a diffusion 
process on a certain family of classical real time paths. By systematically incorporating the phase fluctuations we 
managed to derive for the first time within a real time semiclassical formalism coherent tunneling dynamics. 

Regarding incoherent decay in the deep tunneling regime, the theory revealed the connection between thermal 
equilibrium and the tunneling rate upon which thermodynamic rate formulas are based. In the semiclassical limit flux 
across the barrier and equilibrium are linked via the intimate relation between Euclidian and real time paths in the 
complex plane mechanics. However, the theory not only reproduces the results of various other semiclassical thermal 
rate theories for higher temperatures, but covers with no further assumptions also the low temperature regime where 
so far other real time methods failed. 

Several questions could not be analyzed in detail in this article: There is first the temperature range around T c /2 
where approximately high and low temperature semiclassics match; we briefly sketched corresponding improvements. 
Second, we only touched in passing the explicit real time dynamics in the transient time domain where for incoherent 
processes the relaxation to a stationary flux occurs. Third, other initial preparations, e.g. to gain correlation functions, 
were out of the scope of this paper. 

Moreover, with the appropriate formalism at hand further extentions are possible. While the dynamics of dissipative 
systems was already studied in the high to moderate temperature range |2l| , p6|] , the low temperature tunneling regime 
is now in principle open for investigations. Of course, a crucial point for all further applications is to develop an 
appropriate numerical algorithm to mimic the "diffusion" in the complex plane. We hope to make progress in this 
direction in the near future. 
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FIG. 1. Loop of stationary imaginary and real time paths in the complex time plane z = u + ia. 

FIG. 2. Semilcassical paths (dashed lines) in the complex plane near the parabolic barrier top. Shaded area contains relevant 
intermediate coordinates qi,q[ reached by fluctuations along arrows. 

FIG. 3. Real time paths in the double well potential with wells at ±q a (dots) for various qj and T = (thin lines). The 
thick line shows a typical fluctuation connecting orbits with different <?/. 

FIG. 4. Phase spaces orbits of complex real time paths in the double well potential at T — 0. In the left picture the real 
part of the orbits is shown for trajectories starting with qf > near the barrier top; the dot indicates the well at q a . In the 
right picture the corresponding imaginary parts are depicted. 

FIG. 5. Diffusion of the crossing point g 1 -™' of a fluctuation path along the real axis (thick lines) for various cases discussed 
in the text. Dots indicate the wells at ±q a that are branch points for the momenta, thin vertical lines the end-coordinates at 
±<7/. Solid lines refer to the forward, dotted ones to the backward propagator, a crossing of a dot a TP. 

FIG. 6. Real time paths in the complex plane for the Eckart barrier. Solid lines show orbits for T < T c , dotted lines orbits 
for T > T c . 
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FIG. 7. Phase space orbits of complex real time paths in the Eckart barrier potential at T <C T c . In the left picture the real 
part of the orbits is shown for trajectories starting with 5/ > near the barrier top. In the right picture the corresponding 
imaginary parts are depicted; the dotted line separates the strips n = and n — 1. 

FIG. 8. Transmission factor P as a function of inverse temperature for an Eckart barrier with a = 12, a = 2nVo/h,uJb- P 
is defined by P = r/r c i with T c i the classical rate, i.e. the high temperature limit to Eq. ([u]). The solid line is the exact 
result. In the left picture the dotted line shows the parabolic result Eq. (|3l|), the dashed line represents Eq. (|7^), and the arrow 
indicates the inverse temperature corresponding to T c . In the right picture the dashed line depicts the result Eq. (|3|) and the 
arrow refers to T c /2. 
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Figure 1 , Ankerhold et al 
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Figure 2, Ankerhold et al 
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Figure 3, Ankerhold et al. 
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Figure 4, Ankerhold et al. 
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Figure 5, Ankerhold et al. 
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Figure 6, Ankerhold et al. 
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Figure 7, Ankerhold et al. 
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Figure 8, Ankerhold et al. 



29 



